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ABSTRACT 

Measurements of the cosmological density parameter using techniques that 
exploit the gravity-induced motions of galaxies constrain, in linear perturbation theory, 
the degenerate parameter combination /? = Q^/b, where the linear bias parameter 
b is the ratio of the fluctuation amplitudes of the galaxy and mass distributions. 
However, the relation between the galaxy and mass density fields depends on the 
complex physics of galaxy formation, and it can in general be non-linear, stochastic, 
and perhaps non-local. The one-parameter linear bias model is almost certainly 
oversimplified, which leads to the obvious question: What is the quantity P that is 
actually measured by different techniques? To address this question, we estimate /3 
from galaxy distributions that are constructed by applying a variety of locally biased 
galaxy formation models to cosmological N-body simulations. We compare the values 
of P estimated using three different techniques: a density-density comparison similar 
to the POTENT analysis, a velocity-velocity comparison similar to the VELMOD 
analysis, and an anisotropy analysis of the redshift-space power spectrum. In most 
cases, we find that (3 estimated using all three methods is similar to the asymptotic 
value of /bfj{R) at large i?, where bfj{R) is the ratio of rms galaxy fluctuations 
to rms mass fluctuations on scale R. Thus, something close to the conventional 
interpretation of /? continues to hold even for complex bias models. Moreover, we find 
that P estimates made using these three methods should, in principle, agree with each 
other. It is thus unlikely that non-linear or scale-dependent bias is responsible for the 
discrepancies that exist among current measurements of P from different techniques. 

Subject headings: cosmology: theory, galaxies: distances and redshifts, methods: 
numerical 
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1. Introduction 

In cosmological linear perturbation theory, there is a simple relation between the peculiar 
velocity and the mass density fields, 

v(x) = — — / 5m(x )— —^d X , (1) 



47r J ""'^"'|x' 

where (5m (x) = p{x.)/p — 1 is the mass density contrast, f{^m) ~ ^^'i is the ratio of the 



average mass density of the universe to the critical density ( [Peebles 1980| ). The differential form 
of this relation is 

V • v(x) = -f{Q^)6m{^) . (2) 

These relations suggest that if we map the mass density and velocity fields, a comparison of the 
two will yield a measurement of Vtm- Unfortunately, we cannot observe the mass distribution 
directly. We can only detect luminous matter, of which galaxies are the basic unit. The relation 
between the galaxy and mass distributions, usually referred to as "bias", depends on the details 
of the galaxy formation process. A complete theory of galaxy formation should predict what 
environments galaxies form in and how they are distributed in space with respect to the mass. 
However, while galaxy formation is one of the most actively pursued fields of theoretical cosmology, 
our current understanding of it is far from perfect. The relation between galaxies and mass is 
often parameterized by the linear bias model 6g = b5m, where 6g is the galaxy density contrast 
and b is the linear bias factor. Equation (^) then becomes 

V-v(x) = -/?<53(x) , (3) 

where /? = f{0,rn)/b. 

Methods that use equation (|3|), or some form of it, to infer the cosmic mass density can only 
directly measure /3, a degenerate combination of 0^. and b. In practice, /3 is usually estimated from 
a comparison between the galaxy density field, inferred from galaxy redshift maps, and galaxy 
peculiar velocities, which require distance measurements to individual galaxies. This estimate is 
made either through a density-density comparison, where the mass density field is predicted from 
the smoothed peculiar velocity field and is then compared to the observed galaxy density field, or 
through a velocity-velocity comparison, where galaxy peculiar velocities are predicted from the 
galaxy density field and are then compared individually to observed galaxy peculiar velocities. 
Alternatively, /3 can be measured solely from galaxy redshift maps by analyzing the anisotropy 
of galaxy clustering produced by redshift-space distortions, using either the power spectrum or 
correlation function in redshift space. Based on these three main approaches, there are many 



techniques that have been developed and used to estimate /? (see [Strauss Sz Willick 1995| for a 
review of these techniques). 

Since the bias relation depends on the complex process of galaxy formation, it is probably 
more complicated than the one-parameter linear bias model. It can, in general, be non-linear, 
stochastic, and possibly non-local ( Dekel fc Lahav 1999[ ). Therefore, it is not obvious exactly 



what information a measurement of f3 contains. Moreover, it is not clear whether the different 
techniques for estimating /? are measuring the same quantity. The primary purpose of this paper 
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is to examine these issues and determine what can be learned about bias from measurements of f). 
We create simulated galaxy distributions with various non-trivial bias prescriptions and use them 
to measure (5 with different existing techniques. The general approach is similar to the one we 
followed in Narayanan, Berlind fc Weinberg (2000| ) (hereafter NBW), where we investigated the 



sensitivity of large-scale structure statistics to bias. 

An additional motivation for this study is that there are large discrepancies in current (3 
measurements made with different methods. For example, analyses of the IRAS 1.2-Jy galaxy 
redshift su rvey ([Fisher et al. 1995|) have yielded fSjjiAS = 0.89 it 0.12 by the density-density 
approach ( ^igad et al. 1998|) , Piras = 0-50 ± 0.04 by the velocity- velocity approach ( Willick fc 



[Strauss 1998| ), and (3iras = 0.52 ± 0.13 by the redshift -space distortion approach (Cole, Fisher & 
Weinberg 1995). Table 1 summarizes current f3 estimates for IRAS-selected samples of galaxies 
made using different techniques. The estimates of /? span the range 0.4 — 1.0, and the most 
extreme measurements differ from each other at the ~ 4(T level. Some authors have suggested that 
these discrepancies could be caused by complexity in the bias relation affecting different methods 
in different ways ( [Wilhck et al. 1997[ ; [Sigad et al. 1998| ; [Dekel fc Lahav 1999| ). For example. 



each technique effectively probes the bias on a different physical scale, so it is plausible that a 
scale-dependent bias relation could cause such discrepancies. In this paper we examine whether 
the variation in measured values of (5 can be naturally explained by complexity of bias. 

We employ a variety of simple, non-linear, local bias prescriptions to create galaxy 
distributions from the outputs of N-body simulations. We then measure P using somewhat 
idealized versions of three different methods: a density-density comparison, a velocity-velocity 
comparison, and an analysis of the anisotropy of the redshift-space power spectrum. We assume 
that galaxy positions and peculiar velocities are known perfectly, and we therefore do not attempt 
to model the random and systematic errors that exist in real data. We are thus able to isolate 
the effects of different types of local biasing on measurements of f3 without worrying about how 
these measurements will be affected by observational errors. The systematic errors that exist in 
real observations are, of course, very important, but they are best dealt with in the context of 
specific data sets. We compare our P measurements to two well-defined functions. The first is 
Pct{R) = /ba{R)- Here, ba-{R) is the bias function defined as ba{R) = ag{R)/arn{R), where 
(Tg{R) and am{R) are the rms fluctuations of the galaxy and mass density fields, smoothed with 
a top-hat filter of radius R. The second function is I3p{k) = f{0,rn)/bp{k), where bp{k) is the 



bias function in Fourier space, defined as bp{k) = y Pg{k)/ Pmik), where Pg{k) and Pm{k) are the 
power spectra of the galaxy and mass distributions. In NBW we demonstrated that, for local bias 
models, ba{R) and bp{k) become scale-independent at large scales. We denote the asymptotic 
large-scale values of PaiR) and Pp{k) as and /3p, respectively. 



2. Models 

We have carried out N-body simulations of three different cosmological models, all 
based on infiation and cold dark matter (CDM). The first is an Vtm = 1, h = 0.5 model 
{h = Hq/100 km s~^ Mpc~^), with a tilted power spectrum of density fiuctuations designed 
to satisfy both COBE and cluster normalization constraints. The cluster constraint requires 



- 4 - 



crsm^m^ ~ 0.55 ( [White, Efstathiou k Frenk 1993D , where asm is the rms amphtude of hnear mass 
density fluctuations in top-hat spheres of radius 8/i~^Mpc. Matching the COBE-DMR constraint 
and (Tgm = 0.55 with h = 0.5 requires an inflationary spectral index n = 0.803 if one incorporates 
the standard inflationary prediction for gravitational wave contributions to the COBE anisotropics 
(see Cole et al. 1997 and references therein). The other two models have = 0.2 and 0.4, with a 
power spectrum shape parameter F = 0.25 (in the parameterization of Efstathiou, Bond & White 
1992 ) and a cluster-normalized fluctuation amplitude agm = 0.55i7~'''^. These two models are 
open models with no cosmological constant (Oa = 0). Since Oa has a negligible effect on peculiar 
velocities at fixed $7.^, our results for the two open models should also hold for flat-l^A cosmologies 
having the same values of and the same matter power spectrum. All simulations were run 
with a particle-mesh (PM) N-body code written by C. Park, which is described and tested by 
Park (19"90| ). Each simulation uses a 400^ force mesh to follow the gravitational evolution of 200^ 
particles in a periodic cube 400/i~^Mpc on a side, starting at 2; = 23 and advancing to z = in 46 
steps of equal expansion factor o. We have run four independent realizations of each of the above 
cosmological models, and all of the results shown in this paper have been averaged over these four 
realizations. 

The bias between galaxies and mass should ultimately be a prediction of a theory of galaxy 
formation. There are currently three theoretical approaches to predicting how galaxies are 
distributed in space with respect to the mass. Semi-analytic models of galaxy formation identify 
virialized dark matter halos in moderate resolution N-body simulations and then populate them 



with galaxies using analytic prescriptions (e.g., Kauffmann et al. 1999; Benson et al. 2000); high 
resolution N-body simulations resolve individual galaxy-sized dark matter halos within larger halos 
(e.g., Kravtsov fc Klypin 199S ; Colin et al. 1999 ; Moore et al. 19991 ); hydrodynamic simulations 
follow the evolution of both dark matter and baryons, including the effects of gas cooling, star 
formation, and feedback, and identify galaxies based on the location of cold baryonic lumps (e.g., 
Cen &: Ostriker 1992; Katz, Hernquist & Weinberg 1992; Katz, Weinberg fc Hernquist 199^ ; 
Blanton et al. 199£ ; Pearce et al. 1999| ). All three approaches have made much progress in the 
last few years, but we are still far from a complete and compelling theory of galaxy formation, 
and even the most ambitious applications of these techniques have explored only a small selection 
of cosmological models and relatively small simulation volumes. Rather than adopt a specific, 
detailed theory of bias, we employ a simple approach to modeling the "generic" effects of bias on 
large-scale structure measurements. We apply plausible biasing schemes to select galaxy particles 
from the mass distributions and check which measurements of f3 are sensitive to the details of 



the biasing scheme and which are robust, extending the approach used by NBW (see also Mann 



Peacock & Heavens 199^ 



We create galaxy distributions by applying various local biasing prescriptions to the mass 
distributions. In these bias prescriptions, the probability of a given mass particle being selected 
as a galaxy depends on the properties (density, geometry, or velocity dispersion) of the mass 
distribution, averaged in top-hat spheres of radius 4/i~^Mpc, centered on that particle. These 
prescriptions are also used and described in more detail by NBW and [Narayanan et al. (1999| ) 
brief, the bias prescriptions are: 



In 



(1) Semi-analytic: An empirical bias prescription derived by Narayanan et al. (1999| ) which 
characterizes the relation between the galaxy and mass density fields in the semi-analytic galaxy 
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formation models of pBenson et al. (2000 ). 



(2) Sqrt-Exp (Square-root Exponential): A bias prescription in which (1 + 6g) on 

+ 6m)e°'^^^^"^\ Here, {l + 6g) and {l + 5m) are the galaxy and mass overdensities, respectively. 
This prescription is intended to be monotonic (for a > 0) and non-linear, but it can be tuned to 
allow galaxies to be either biased or anti-biased with respect to the mass, and it yields a non-trivial 
bias relation even when = 1. 

(3) Power-law: A bias prescription in which (1 -|- 6g) oc (1 -|- (5m)"- This is similar to the bias 
relation suggested by Cen fc Ostriker (1993| ), based on hydrodynamic simulations. 



(4) Threshold: A bias prescription in which galaxies do not form below some mass density 
threshold, and they form with equal efficiency at densities above that threshold. 

(5) Sigma: A bias prescription in which galaxies are selected with equal probability from mass 
particles that have a velocity dispersion av greater than some threshold value. 

(6) Sheet: A bias prescription that selects galaxies in regions where the mass distribution is 
planar. We compute the eigenvalues (A3 > A2 > Ai) of the moment of inertia tensor in spheres 
of radius 4/i~^Mpc around each mass particle and select galaxies to be those mass particles that 
have the highest ratio of A3/A1. 

(7) High-z: A bias prescription identical to the Power-law bias, except galaxies are selected 
based on the mass distribution at redshift 2 = 3. This prescription is intended to produce a 
biasing relation that is stochastic on the 4/i~^Mpc scale at 2; = 0, though still biased in the mean. 

Each of the bias models contains one tunable parameter that is adjusted so that the resulting 
galaxy distribution has an rms fluctuation, in top-hat spheres of radius 12/i~^Mpc, of CJ12 ~ 0.7. 
The values of &(t(12) for the galaxy distributions are 1.7, 1.0, and 0.67 for VLm = 1-0, 0.4, 
and 0.2, respectively. In all the analyses in this paper, we form the mass and galaxy density 
fields by cloud-in-cell (CIC) binning the particle distributions onto a 200'^ grid. We create a 
volume-weighted, smoothed velocity field from the discrete galaxy peculiar velocities using the 



method of Babul et al. (1994 ). Specifically, we first form the momentum field by CIC-binning 
the momentum of every galaxy onto a 200'^ grid. We smooth this momentum field with a 
Gaussian filter of radius Ri = Rs/2 and divide it by a similarly smoothed density field to form a 
mass-weighted smoothed velocity field. We then smooth this velocity field with another Gaussian 
filter of radius R2 = {R?s ~ -R?)^^^) so that the effective smoothing radius is Rg- Because the 
second smoothing dominates over the first, the final velocity field is volume- weighted rather than 
mass- weighted. 



3. Density-Density Comparison 

Density-density comparisons make direct use of equation (^) to measure /3. The whole process, 
however, involves many steps. Radial peculiar velocities must first be computed for a sample of 
galaxies that have both redshift and (redshift-independent) distance measurements. Due to the 
large uncertainties present in galaxy distance measurements, individual galaxy radial peculiar 
velocities are poorly known. The resulting radial velocity field must therefore be smoothed over 
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a large scale to reduce the velocity errors. The full 3-dimensional velocity field is then inferred 



from the radial velocity field using the POTENT method proposed by [Bertschinger k. Dekel 



19891) , under the assumption that galaxy velocities trace the gravitational potential field of the 



underlying mass distribution. Finally, the divergence of this smoothed 3-dimensional velocity field 
yields the linear-theory prediction for the smoothed mass density field. On the other side of the 
density-density comparison, the real space galaxy density field is obtained from a galaxy redshift 
catalog, after correcting the galaxy redshifts for peculiar velocities estimated using the redshift 
space density field itself. The real-space galaxy density field is then smoothed on the same scale 
as the velocity field. After all of these steps are taken, the observed galaxy density field may be 
directly compared to the predicted mass density field. As shown by equation (P), the slope of this 
comparison yields fi. 



In practice, a non- linear generalization of equation (y) is used to determine (3 (e.g., Dekel 



et al. 1993|; Hudson et al. 1995| ; |Sigad et al. 1998| ; pekel et al. 19991 ). For example, ^igad et al 



(1998 ) measure (3 from a comparison of the IRAS 1.2-Jy galaxy density field, 5iras^ to the mass 
density field derived from Mark III peculiar velocities. They predict the mass density field using 
the non-linear approximation 

6poT = -(1 + ei)/-iV • V + (1 + e2)/-'A2 + (1 + e3)/-'A3 , (4) 

where A2 and A3 are second and third order terms that involve sums of double and triple 
products of partial derivatives of the velocity field, and the three coefficients ei, £2 and £3 are 
empirically determined from a family of CDM N-body simulations. In addition, Sigad et al. (199^ ) 



use a non-linear approximation to convert the IRAS galaxy redshift map to real space. Since 
equation (^) contains higher powers of f(Qrn), it is not possible to measure /3 simply by fitting a 
line through the Sjras vs. 5pot relation, ^igad et al. (1998|) first assume a value for $7^ and then 
fit a line through the relation to get Their estimate for [3 is then equal to f{^rn)/b. Since the 
initial choice of 0^ does not seem to affect the final /? measurement, it appears that it is indeed a 
degenerate combination of and b, even in this mildly non-linear regime. 

Aside from addressing non- linear effects, density-density analyses must deal with a host 
of problems that arise due to the imperfect nature of observational data sets. In particular, 
the peculiar velocity errors, when smoothed, give rise to inhomogeneous Malmquist bias, which 
is difficult to correct for. Moreover, density-density analyses must deal with issues such as 
non-uniform sampling and imperfect survey boundaries in both the galaxy redshift, and peculiar 
velocity data. Sparse sampling of the velocity field, for example, makes it particularly hard to 
smooth the velocity data in an unbiased way. Much work has been done recently to understand 
and control all of these problems ( pekel et al. 199^ ). Current density-density comparisons smooth 
the velocity field with a Gaussian filter of radius 12/i~^Mpc and yield rather high values of (3 
compared to other techniques (see Table 1 for a summary of current measurements) . The analysis 



of ISigad et al. (1998|) yields f3iRAs = 0.89 ± 0.12. 



We wish to focus directly on the infiuence of bias on density-density comparisons, and we 
therefore adopt a simpler approach to analyzing our numerical data sets. We assume perfect 
knowledge of both the real-space positions and the velocities of galaxies in our simulation volume. 
We then smooth the galaxy velocity field with a Gaussian filter of some radius R and compare 
its divergence to the galaxy density field smoothed at the same scale. We fit a line of form 
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— (V • v) = f36g + C to the density-density relation in the region —0.5 < Sg < 0.5 in order to 
estimate /3. Figures 1 and 2 show this procedure when the velocity and density fields are smoothed 
at 12/i~^Mpc for the = 1-0 and 0.2 models. In each panel we show the mean relation between 

— (V • v) and 6g for a given bias prescription, along with its la scatter, and the best-fit line to 
this relation (we have omitted the Sigma and High-z bias prescriptions because they produce 
results that are nearly identical to the Threshold and Power-law prescriptions, respectively). The 
slope /3fit of the best-fit line is indicated at the lower right corner of each panel. Also shown, for 
purposes of comparison, is the asymptotic large-scale value of Pa{R). Table 2 summarizes the 
results for all the biasing models (including the Sigma and High-z bias models) for each of the 
three cosmological models. 

For = 1.0 and 0.4, all the bias models, with the exception of the Sqrt-Exp and Sheet 
models, have best-fit estimates of /? that are in excellent agreement with the asymptotic large-scale 
value Pa, deviating not more than 4% from this value. The Sqrt-Exp model yields a best-fit value 
that is substantially higher than whereas the Sheet model yields a /3 estimate that is somewhat 
lower than /?o-. However, the strong degree of non-linearity in the — (V • v) vs. 5g relation for 
both of these models is striking in appearance and can probably be ruled out with existing data 



(e.g., Sigad et al. 199S|). For = 0.2, the situation is slightly different. The Semi-analytic 



and Power-law models yield P estimates that underestimate Per by ~ 10%, whereas the Sqrt-Exp 
model yields a P estimate that overestimates Pa by 8%. Of course, analyses of real data sets must 
deal with the additional complications of non-linearity and statistical biases, but these have been 
addressed in the papers cited above and depend in detail on properties of the data sets themselves. 
Our results show that once these challenges are met, the value of P derived from POTENT-like 
analyses should be close to the asymptotic value of Pa, for a fairly broad range of assumptions 
about the form of biasing. 



4. Velocity- Velocity Comparison 

Velocity-velocity comparisons make use of equation (|l]) to measure p. Roughly speaking, the 
observed galaxy redshift distribution is used to predict the peculiar velocities of individual galaxies 
using equation (||), with /(J7m) replaced by an assumed value for p. These velocity predictions 
are then compared directly to peculiar velocity measurements made using redshift-independent 
distance measurements. The best estimate of P is that for which the predicted and measured 
galaxy velocities show the best agreement. 

Two distinct velocity-velocity methods have been used in recent measurements of p. The 
first method, described by Willick et al. (1997| ), requires that the galaxy distribution be converted 



from redshift space to real space before it is used to predict peculiar velocities. This is a slightly 
tricky step, since the peculiar velocities are both the final product and a required intermediate 
ingredient in this process. Once the galaxy positions are corrected for redshift-space distortions, 
the galaxy density field must be smoothed in order to suppress the effects of non-linear evolution 
and shot noise. Equation (|l]) is then used to compute the predicted velocity field from the 
smoothed galaxy density field, assuming a value for /?. Finally, a maximum likelihood analysis 
(VELMOD) is used to find the value of P that produces the best match between the predicted and 
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observed galaxy velocities, assuming that an individual galaxy's velocity is the sum of the linear 
theory prediction and an uncorrelated "thermal" velocity. In practice, there are many technical 
difficulties involved in this process. To compute the predicted velocity field, the galaxy density 
field must be integrated over all space, as shown by equation (|l|). Consequently, any systematic 
problems in the galaxy redshift data (such as empty regions) will affect the predicted velocity field 
everywhere. Moreover, the true velocity at any location will be affected to some extent by density 
features that are outside the volume sampled by a redshift survey. [Wilhck et al. (19971) address 



this issue by adding to the predicted velocity field a quadrupole term that models the tidal field 
arising from density features external to the volume probed by the IRAS 1.2-Jy redshift survey. 
The current VELMOD analysis compares the velocities predicted from the galaxy density field 
of the IRAS 1.2-Jy redshift survey, smoothed with a Gaussian filter of radius 3/i~^Mpc, to the 
peculiar velocities in the Mark III catalog and yields Pibas = 0.5 ± 0.04 ( Willick fc Strauss 1998D . 



The second velocity- velocity method, described by [Nusser &: Davis (1994 ) (hereafter ND94) 



and first implemented by Davis, Nusser & Willick (1996), compares the galaxy density field to the 
observed peculiar velocity field directly in redshift space. Another important difference between 



this method and VELMOD is that Davis et al. (1996 ) do not compare predicted and observed 



velocities of individual galaxies. Instead, they expand both the predicted velocity field that is 
derived from the redshift-space galaxy density field and the velocity field that is measured from 
redshift-independent distances into a set of orthogonal modes. They then perform a mode-by-mode 
comparison of the two fields and determine the value of f3 for which the best match is obtained. 
The ND94 method has been applied to estimate [iiRAS by combining the 1.2-Jy redshift survey 
with the Mark III peculiar velocity sample, yielding Piras = 0.6 = fc 0.2 (pavis et al. 19961 ), with a 
surface-brightness fluctuation sample, yielding (Siras = ^■^'^-o oe ( |Blakeslee et al. 1999|) , and with 



a type la supernova sample, yielding Piras = 0.4 it 0.15 ( [Riess et al. 1997 ). Most recently, [Nusser 



et al. (20001 ) have applied this method to the Point Source Catalog redshift survey (PSCz) and the 



ENEAR peculiar velocity sample to obtain Piras = 0.5 ± 0.1. 

In order to investigate the effects of bias on velocity-velocity measurements of (3, we perform 
a simplified VELMOD-like analysis. As in the density-density comparisons, we assume perfect 
knowledge of both the real-space positions and velocities of galaxies in our simulation volumes. 
We smooth the galaxy density field with a Gaussian filter of radius R and use equation (|l|) to 
predict the velocity field. We then interpolate to galaxy positions to find predicted velocities 
(vpred) for all the galaxies. Finally, we compare these with the true galaxy velocities (vtrue)- The 
slope of the best-fit line to the vtrue vs. Vpred relation is our measurement of p. We do not include 
a constant offset term in our fit as we did for the density-density analysis, but it would make 
little difference to our results because the vtrue vs. Vpred relations are very close to linear. The 
offset term is important in the density-density analysis because of the non-linearity of the Jtrue vs. 
(5pred relations. Figures 3 and 4 show this procedure when the galaxy density field is smoothed 
with a Gaussian filter of radius 3/i~^Mpc for our = 1-0 and 0.2 models. In each panel we 
show the mean relation between vtrue and Vp^ed for a given bias prescription, along with its la 
scatter, and the best-fit line to this relation (as before, we have omitted the Sigma and High-z 
bias prescriptions because they produce results that are nearly identical to the Threshold and 
Power-law prescriptions, respectively). The slope Pat of the best-fit line is indicated at the lower 
right corner of each panel. Also shown, for purposes of comparison, is the asymptotic large-scale 
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value of (3(j{R). Table 2 summarizes the results for all the biasing models (including the Sigma 
and High-z bias models) for each of the three cosmological models. 

For all three cosmological models, all the bias models, with the exception of the Jim = 1-0 
Sheet model, have best-fit estimates of (3 that agree well with the asymptotic large-scale value fi^, 
deviating not more than 11% from this value. The Sheet model yields a /3 estimate that is 33% 
lower than (3^- 

For smoothing lengths larger than 5/i~^Mpc, we discovered a systematic error that affects (3 
estimates made using the VELMOD method, even in the case where galaxies trace mass. It arises 
because a smoothed quantity (vpi-ed) is compared to an unsmoothed quantity (vtrue)) and the errors 
in the predicted velocities are thus correlated with the predicted velocities themselves. This issue 
is fully explored in a previous paper ( [Berlind, Narayanan &: Weinberg 200[1| ). The mode- by- mode 
method of ND94 should not be affected by this systematic error because it effectively smooths 
both the velocity and the density fields in the same manner. For a smoothing length of 3/i~^Mpc, 
the infiuence on VELMOD analyses is small, at least for the cosmological models considered here. 



5. Anisotropy Of Redshift- Space Clustering 

Both density-density and velocity- velocity methods require measurements of galaxy peculiar 
velocities in addition to redshifts. Peculiar velocity measurements always contain large statistical 
and systematic uncertainties because they rely on distance-indicator relations that have substantial 
intrinsic scatter as well as uncertain calibrations and environmental dependencies. Moreover, 
these methods usually require the galaxy density field in real space and therefore rely on a 
conversion of the observed redshift-space galaxy density field. There is an entirely different 
approach to estimating (5 that does not suffer from either of these problems, though it does 
require very large redshift samples for effective application. It takes advantage of the fact that 
line-of-sight distortions in redshift space are caused by the same velocities that density-density 
and velocity-velocity methods must measure independently. It is possible to estimate /3 simply by 
analyzing these redshift-space distortions, or, more specifically, by measuring the anisotropy of the 
galaxy clustering in redshift space. 
I 1 

Kaiser (1987) showed that, in the linear regime, the redshift-space galaxy power spectrum 
[Pg{k,^)] is related to the real-space galaxy power spectrum [Pj^^k)] by 

p^{K^,) = {l + p^J?fp^{k), (5) 

where is the cosine of the angle between the line-of-sight and the wave- vector k of a fiuctuation 
in the galaxy density field. This equation is derived using the plane-parallel approximation, i.e, 
in the case where the volume probed is distant enough that all lines of sight to it are effectively 
parallel. Equation (||) reveals that the redshift-space power spectrum is anisotropic. The power 
spectrum of fluctuations along the line-of-sight is amplified by the amount (1 + /3)^ with respect 
to the power spectrum in real space, whereas the power spectrum of fluctuations perpendicular to 
the line-of-sight is not affected at all. The redshift-space power spectrum (^) has been shown to 
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reduce to a sum of monopole (/ = 0), quadrupole (1 = 2), and hexadecapole {I = 4) terms, 

if (fc, /i) = Pi{k)Lo{ij) + Pi{k)L2{fi) + Pf (fc)L4(/i), (6) 

where Li{fj,) are Legendre polynomials. The first two moments are related to the real space galaxy 
power spectrum by 

Piik) = + (7) 

(Cole, Fisher & Weinberg 1994). Hamilton (1992) derived a similar set of relations for the 



multipole moments of the redshift-space galaxy correlation function. Equations @ and (0) provide 
us with two different ways to measure /? from the redshift space power spectrum P^{k,^) (from 
now on we will drop the subscript g in the galaxy power spectrum); using P^ (k)/ P^{k), the ratio 
of the angle-averaged redshift-space power spectrum (monopole) to the real-space power spectrum, 
or using P2{k)/P(){k), the ratio of the quadrupole and monopole moments of the redshift-space 
power spectrum. These ratios are, in principle, measurable and in linear theory they are related 
to (3 by 

^^^^^ - 1 + + (8) 



pR{k) 3"^ 5' 

^^^^^^ ' (9) 



Poik) (1+2/3+1/32) 



These ratios yield an estimate of /? at each wavenumber k. However, non-linearity in 
the velocity and density fields produces distortions in an opposite sense to the linear theory 
predictions, leading to systematically lower estimates of /3 even on scales as large as 50/i~^Mpc. 
Hence, it is essential to model the non-linearities accurately in order to estimate /? using currently 



available redshift surveys. Cole et al. (1995 ) estimated (3 from P2{k)/PQ(k) by assuming an 
exponential velocity distribution model, in which galaxies have uncorrelated small scale peculiar 
velocities drawn from an exponential distribution in addition to their linear theory velocities. 



Applying this method to the IRAS 1.2-Jy redshift survey, they found Piras = 0.52 it 0.13. Fisher 



&; Nusser (1996 ) modeled the non-linearity in P2{k) / PQ{k) by using the Zel'dovich approximation 



( Zel'dovich 1970| ), thus assuming that the scale dependence in this ratio is caused by coherent. 



rather than random, non-linear motions. They found Piras = 0.6 it 0.2. Finally, [Hatton &: Cole 



19991) proposed and tested an empirical model for the non-linearity in P2{k) / Po{k) by examining 



the scale dependence of this ratio in a large number of N-body simulations spanning a broad range 
of cosmological parameters. This model is more general than the previous ones because it is based 
on fully non-linear N-body data. 

Equations (|8|) and only hold in the linear regime and the plane-parallel approximation. In 
order to use these relations to measure (3, we must measure P^{k,^) in volumes that are both 
large (so that they contain fluctuations large enough to be in the linear regime) and far away 
(so that all lines of sight to a single volume are approximately parallel). These constraints make 
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it difficult to accurately measure /3 from the redshift surveys that exist today. Ongoing surveys, 
such as the 2dF redshift survey and Sloan Digital Sky Survey (SDSS), will be much better suited 
to this purpose because they will probe very large volumes, although even then non-linearity will 
be important. To avoid the need for the distant observer approximation, Fisher, Scharf & Lahav 
(1994) measured /3 from the IRAS 1.2-Jy redshift survey by expanding the galaxy redshift-space 
density field into spherical harmonics and maximizing the likelihood that an assumed real-space 
galaxy power spectrum would yield that specific set of harmonics. The free parameters in this 
spherical harmonic analysis (SHA) are /3 and T, the latter determining the shape of the real-space 
power spectrum. Fisher et al. (199^ obtained Piras = 0-94 ± 0.17. More recently, [Ballinger et 



al. (200q) used the SHA method to measure 

PiRAS = 0.40 ± 0.10 for the PSCz, and [Hamilton 



Tegmark &: Padmanabhan (200d| ) developed an optimal way to apply the SHA method to the 



same survey and measured Pjras = 0-4llo'i2- Table 1 gives a summary of (3 measurements from 
redshift-space distortions in IRAS selected galaxy redshift surveys (also see [Hamilton 199^ for a 
review of such measurements) . 

For our simulation analyses, we estimate (3 from the anisotropy of the redshift-space power 
spectrum using an approach that is similar to that of Cole et al. (1995| ), but somewhat simplified. 



We take the line-of-sight direction to be a Cartesian axis of the simulation cube, implicitly 
assuming that the whole simulation volume is far enough away to satisfy the distant observer 
approximation. We measure P^{k,fi) and P^{k) using a Fast Fourier Transform (FFT), and 
we extract the multipole moments by fitting equation (^) to P^{k,fj,). We compute the average 
values and la uncertainties of the ratios P^{k)/P^{k) and P2{k)/PQ{k) from the four independent 
realizations of each galaxy distribution. The points in Figures 5 and 6 represent the functions 
/5(A) obtained from solving for (5 in equations (^)and @, respectively, where A = lir/k. These 
functions asymptote to a constant value only on large scales, since equations ^ only hold under 
the assumption of linear theory. We estimate (3 for each bias model using two methods. (1) We 
use the [Cole et al. (1995| ) exponential velocity distribution model to model the non-linearity in 
P^{k)/P^{k), and we estimate /3 by fitting this model to our measured P^{k)/P^{k). We perform 
the fit only using modes that correspond to scales A > 20/i~^Mpc, and the fitting parameters 
are /Sgt and cr^ (the small-scale velocity dispersion). The thick lines in Figure 5 represent the 
resulting fits for = 1-0. (2) We use the [Hatton fc Cole (1999[ ) empirical model to model the 
non-linearity in P2{k)/PQ{k), and we estimate f3 by fitting this model to our measured P2(^)/-Po(^)- 
As before, we perform the fit only using modes that correspond to scales A > 20/i~^Mpc, and 
the fitting parameters are /3fit and fcni (the wavenumber that corresponds to the scale at which 
the ratio P2{k) / P^^k) is equal to zero). The thick lines in Figure 6 represent the resulting fits 
for = 1-0. The thin lines in Figures 5 and 6 show the function (3p{X) = r2^^/6p(A), where 



bp{\) = \J Pg{k) / Pm{k) at k = 27r/A. In both figures, the best-fit estimates /^gt are indicated at 
the lower right corner of each panel. Also shown, for purposes of comparison, are the asymptotic 
large-scale values of (3(y{R). Table 2 summarizes these values and also includes results for the 
Sigma and High-^; bias models, which are omitted from Figures 5 and 6, and for the low fim 
cosmological models. 

For all three cosmological models and all bias models, with the exception of the Vim = 1-0 
Sqrt-Exp and Sheet models, the best-fit estimates of (3 made using the redshift-to-real-space ratio 
P^{k)/P^{k) underestimate /3o- by 10 — 20%. This underestimate occurs even for the unbiased 
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(mass) models, indicating that it arises because the exponential velocity distribution model for 
the non-linear behaviour of P^{k)/P^{k) is not accurate at this level. Cole et al. (1995| ) reached a 
similar conclusion about the accuracy of this non-linear model, in a slightly different manner. The 
patton fc Cole (1999| ) empirical model does a better job fitting the quadrupole-to-monopole ratio 
P2{k)/ Po{k). For Jim = 1.0, all of the bias models, with the exception of Sqrt-Exp and Sheet, 
yield /? estimates that underestimate by 3 — 8%. For low all of the bias models, again with 
the exception of Sqrt-Exp, yield /? estimates that agree well with /Jq-- 



6. Summary 

Figures 7, 8 and 9 summarize all of our results for = 1-0, 0.4, and 0.2, respectively. 
Each panel corresponds to a particular biasing prescription and shows /? estimates, as a function 
of smoothing length, made using a density-density comparison (circles) and a velocity-velocity 
comparison (squares). Also shown, at arbitrary scales, are the /3 estimates made from the 
anisotropy of redshift-space clustering (triangles). Solid points represent the smoothing scales 
that correspond to recent observational estimates of /3: the density-density (POTENT-like) 
estimate made with a 12/i^^Mpc Gaussian smoothing (solid circle; cf. [Sigad et al. 199^ ), the 
velocity-velocity (VELMOD-like) estimate made with a 3/i~^Mpc Gaussian smoothing (solid 
square; cf. Willick et al. 199^) , and the redshift-space anisotropy estimate made using the 



quadrupole-to-monopole ratio of the redshift-space power spectrum, P2{k)/ PQ{k), marked at an 



arbitrary scale (solid triangle; cf. Cole et al. 1995 ). Also shown, for comparison, is the function 
Pa{R) (solid line). 

These results allow us to answer the question "what is /5?". In most cases, to a fairly good 
approximation, the quantity f3 estimated by these methods is close to the ratio Q^/ba, where bo- is 
the asymptotic, large-scale value of the function ag{R) / am{R)- This is true for all three methods 
and for all three cosmological models. This reassuringly simple answer is just what one might have 
hoped for. In our models, the bias relation in non-linear and bias factors like b(j{R) and bp{k) 
are scale-dependent on small scales. Nonetheless, the value of (3 has a well-defined, intuitively 
sensible meaning that holds for a wide range of such models. Furthermore, measurements of [3 
from density-density comparisons, velocity-velocity comparisons, and analyses of the anisotropy of 
redshift-space clustering, should, in principle, yield consistent results. 

There are some exceptions to this simple characterization. The clearest is the Sqrt-Exp bias 
model, which has a strongly non-linear relation between galaxy and mass density. This model 
exhibits the greatest scale dependence of bcj{R) and bp{k), and different methods of estimating /3 
give very different results. However, the strikingly non-linear shape of the — (V • v) vs. 6g relation 
for this model can be ruled out by current POTENT analysis of observational data (|Sigad et al. 



1998), so a bias relation with such pathological effects is probably unrealistic. Another exception 
is the Sheet bias model. Here the different methods yield measurements of (3 that agree well with 
each other, but they do not match the large-scale value Pa- 

Our results imply that non-linear local bias is unlikely to account for the large discrepancies 
that exist between current observational estimates of Piras using different techniques. We 
have also performed these analyses for the non-local bias model considered by NEW, where the 
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relation between galaxy and mass density on the 4/i-~ Mpc scale is modulated by the larger scale 
environment. Although the model does not have a clearly defined asymptotic value of Pa, the 
different /3 estimation methods nonetheless give results that are consistent with each other. We 
therefore conclude that the discrepancies among current observational estimates of /3 probably 
arise from the interaction between the systematic errors in the observational data sets (the redshift 
and peculiar velocity catalogs) and the specific details of the various analysis methods. If the 
observational challenges can be overcome, measurements of /3 by these methods can yield a useful, 
physically meaningful quantity: divided by a bias factor that characterizes the ratio of rms 
galaxy and mass fluctuations, on large scales. 

We note with great sadness the passing of Jeff Willick, who has been a leader in this field for 
more than a decade. We arc grateful to Jeff and to Scott Gaudi and Michael Strauss for helpful 
input and comments. This work was supported by NSF grant AST-9802568. AAB and VKN were 
supported by Presidential Fellowships from the Graduate School of The Ohio State University 
during phases of this project. 
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Table 1. Current (5 measureraents for IRAS selected galaxies 



Data 



Method 



Paper 



/3 



IRAS 1.2-Jy + Mark III 



IRAS 1.2-Jy + SNIa 
IRAS 1.2-Jy + SBF 
IRAS 1.2-Jy 



IRAS QDOT 

PSCz + ENEAR 
PSCz 



POTENT 

VELMOD 

ND94 

ROBUST 

ND94 

ND94 

P2{k)/Po(k) 
P2{k)/Poik) 
SHA 

p^{k)/p^{k) 

P2{k)/Po{k) 
ND94 
SHA 
SHA 



Sigad et al. (1998|) 

Strauss (19981 ) 



Willick 



Davis et al. ( 1996|) 
Rauzy k Hendry (200q) 



Riess et al. 



Blakeslee et al. (1999| ) 



Cole et al. (1995|) 



^T997|) 



Fisher fc Nusser (1996| ) 



Fisher et al. (1994| ) 



Bromley, Warren &: Zurek (1997| ) 
Peacock fc Dodds (19941) 
Cole et al. (19951 ) 



Nusser et al. (2000|) 
Ballinger et al. (2000|) 



Hamilton, Tegmark &: Padmanabhan (20001 ) 



0.89 ±0.12 
0.50 ± 0.04 
0.60 ± 0.20 
0.60 ±0.125 
0.40 ±0.15 



+0.10 
-0.06 



0.42 
0.52 ±0.13 
0.60 ±0.20 
0.94 ±0.17 



+0.40 
-0.30 



0.80 
1.00 ±0.20 
0.54 ± 0.30 
0.50 ±0.10 
0.40 ±0.10 
41+0-13 

'-'■^^-0.12 



Note — We have listed only the most recent estimates for a given data set from each group. 
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Table 2. Estimates of /3 from the biased models, using different techniques 







POTENT 


VELMOD 




-| — V / 1 \ 1 -I— » / 1 \ 

P2{k)/Po{k) 






> 


> 






Q = 1.0 












Mass 


1.00 


0.92 


0.94 


0.80 


0.93 


Semi-analytic 


0.61 


0.61 


0.56 


0.54 


0.57 


Sqrt-Exp. 


0.62 


0.84 


0.69 


0.64 


0.98 


Power-law 


0.60 


0.60 


0.57 


0.52 


0.55 


Threshold 


0.57 


0.55 


0.56 


0.47 


0.55 


Sigma 


0.54 


0.53 


0.56 


0.46 


0.52 


Sheet 


0.66 


0.53 


0.44 


0.46 


0.52 


High-z 


0.60 


0.58 


0.55 


0.52 


0.58 


n = 0.4 












Mass 


0.58 


0.56 


0.54 


0.50 


0.64 


Semi-analytic 


0.57 


0.55 


0.54 


0.49 


0.59 


Sqrt-Exp. 


0.77 


0.93 


0.80 


0.69 


0.76 


n = 0.2 












Mass 


0.38 


0.36 


0.35 


0.32 


0.43 


Semi-analytic 


0.52 


0.46 


0.48 


0.41 


0.51 


Sqrt-Exp. 


0.66 


0.71 


0.64 


0.59 


0.57 


Power-law 


0.54 


0.49 


0.48 


0.45 


0.52 



Note — POTENT results are shown for a 12/i~^Mpc Gaussian smoothing. VELMOD results are 
shown for a 3h~^Mpc Gaussian smoothing. For $7^ = 1.0, all uncertainties are ~ 0.005. For low 
rtjji, uncertainties are ~ 0.002, except for Sqrt-Exp models, where they are ~ 0.04 for = 0.4 

and ~ 0.02 for Urn = 0.2. 
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Fig. 1. — Density-density (POTENT-like) (3 estimates for our $7^ = 1-0 cosmological model. Each 
panel shows the mean measured relation between — (V • v) and 5g (thin solid curve), when the 

velocity and density fields are smoothed with a Gaussian filter of radius 12/i~^Mpc, for a particular 
biasing prescription (we have omitted the Sigma and High-z prescriptions because they produce 
results that are identical to the Threshold and Power-law prescriptions, respectively). The plotted 
relation represents the average over four independent simulations, and error bars indicate the la 
dispersion in this mean relation. The outer two thin solid curves represent the locus of points that, 
for each bin in 5g, enclose 80% of the points in that bin. Also shown is the best-fit line to the 
— (V • v) vs. 5g relation (thick solid line). The slope of this line, denoted by /3fit, is listed in the 
lower right corner of each panel. Also listed for comparison is the asymptotic, large-scale value of 

PAR). 




Fig. 2. — Density-density (POTENT-like) /3 estimates for our Qm = 0.2 cosmological model. Refer 
to Fig. 1 for a complete description of the components of this figure. 
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Fig. 3. — Velocity-velocity (VELMOD-like) P estimates for our = 1.0 cosmological model. 
Each panel shows the mean measured relation between Vtrue and Vpred (thin solid curve), when 

the predicted velocities arc computed from the galaxy density field smoothed with a Gaussian 
filter of radius 3/i^^Mpc, for a particular biasing prescription (we have omitted the Sigma and 
High-z prescriptions because they produce results that are identical to the Threshold and Power- 
law prescriptions, respectively). The plotted relation represents the average over four independent 
simulations, and error bars indicate the la dispersion in this mean relation. The outer two thin 
solid curves represent the locus of points that, for each value of Vprcj, enclose 80% of the galaxies. 
Also shown is the best-fit line to the vtrue vs. Vpred relation (thick solid line). The slope of this 
line, denoted by /3fit, is listed in the lower right corner of each panel. Also listed for comparison is 
the asymptotic, large-scale value of Pa{Fl-)- 




Fig. 4. — Velocity-velocity (VELMOD-like) (3 estimates for our Qm = 0.2 cosmological model. 
Refer to Fig. 3 for a complete description of the components of this figure. 
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(a) Mass (b) Semi — analytic (c) Sqrt-Exp 




A[h-iMpc] A[h-iMpc] A[h-iMpc] 

Fig. 5. — Estimates of /3 from P^{k)/P^{k), the ratio of the redshift- and real-space power spectra, 
for our = 1.0 cosmological model. Each panel shows an estimate of /5 for a particular biasing 
prescription (we have omitted the Sigma and High-z prescriptions because they produce results that 
are identical to the Threshold and Power-law prescriptions, respectively). The points show /9(A) as 
estimated from P^{k)/P^{k) using linear theory and the plane-parallel approximation (eq. P). The 
la errors in the mean /?(A) are computed using the dispersion among four independent simulations, 
divided by \/3. The thin line represents the function /3p(A) = /bp{X), where 6p(A) is the 

bias function, defined as bp{k) = ^ Pg{k)/ Pm{k), where Pg{k) and Pm{k) are the power spectra of 
the galaxy and mass distributions, respectively. The thick line represents a fit of the exponential 
velocity distribution model to /3(A). The fit yields a global estimate of /3fit which is listed at the 
bottom right corner of each panel. Also listed for comparison is the asymptotic, large-scale value 
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(a) Mass (b) Semi-analytic (c) Sqrt-Exp 




U I I L_J tu I I i_J hil I 

10 100 10 100 10 100 



A[h-iMpc] A[h-iMpc] A[h-iMpc] 



Fig. 6. — Like Fig. 5, except that /3(A) is estimated via equation (^) from P2{k)/Po{k), the 
quadrupole-to-monopole ratio of the redshift-space power spectrum, and thick hues represent the 
fit of the patton fc Cole (1999| ) non-hnear model to /3(A). 
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(a) Mass 



(b) Semi-analytic 



(c) Sqrt-Exp 
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(d) Power-law 
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(f) Sigma 



/S(R) 



0.5 




0.5 



(g) Sheet 



(h) High-z 



/S(R) 




R [h-iMpc] 



R [h-iMpc] 



Fig. 7. — Comparison of (3 estimates from different methods, as a function of scale, for our Q^n = 1-0 
cosmological model. Each panel shows l3{R) for a particular biasing prescription. The solid line 
represents the function Prj{R) = /ba{R), with ba-{R) = ag{R) / cjm{R) for a Gaussian filter of 
radius R. Circles represent density-density (POTENT-like) (3 estimates, derived by fitting a line to 
the measured relation between — (V • v) and 5g when the velocity and density fields are smoothed 
with a Gaussian filter of radius R (see Fig. 1). The (3 estimate at a 12/i~^Mpc smoothing radius, 
which corresponds to current POTENT measurements, is highlighted as a solid circle. Squares 
represent velocity-velocity (VELMOD-like) [3 estimates, derived by fitting a line to the measured 
relation between vtrue and Vpred when the predicted velocities Vpred are estimated from the galaxy 
density field smoothed with a Gaussian filter of radius R (see Fig. 3). The (3 estimate at a 3/i~^Mpc 
smoothing radius, which corresponds to current VELMOD measurements, is highlighted as a solid 
square. Triangles show estimates of (3 derived from the anisotropy of the redshift-space power 
spectrum. These estimates are not scale-dependent and are marked at arbitrary values of R. The 
open triangle shows the estimate of [3 derived from fitting the exponential velocity distribution 
model to P^{k)/P^{k) (see Fig. 5). The solid triangle shows the estimate of j3 derived from 
fitting the [Hatton fc Cole (1999| ) non-linear model to P2{k) / Po{k) (see Fig. 6). In all cases, each 
point represents the average over four independent simulations, and the errorbar represents the la 
uncertainty in the mean. 
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Fig. 8. — Comparison of (3 estimates from different methods, as a function of scale, for our = 0.4 
cosmological model. Each panel shows P{R) for a particular biasing prescription. Refer to Fig. 7 
for a complete description of the components of this figure. 
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Fig. 9. — Comparison of /3 estimates from different methods, as a function of scale, for our = 0.2 
cosmological model. Each panel shows P{R) for a particular biasing prescription. Refer to Fig. 7 
for a complete description of the components of this figure. 



